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FOREWORD 

This report documents numerical results obtained by per- 
sonnel of the Computational Mechanics Section at Lockheed's 
Huntsville Engineering Center under Contract NAS8-36284, 
"Turbine Rotor/Stator Flowfield Analysis." 

The NASA-MSFC Contracting Officer's Representative for 
this study is Dr. Helen V. McConnaughey , ED32. 
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ABSTRACT 

A numerical study of the unsteady aerodynamic and thermal environment 
associated with axial turbine stages is presented. Computations are performed 
using a modification of the R0T0R1 rotor/stator interaction code developed by 
M.H. Rai. Two different turbine stages are analyzed: the first stage of the 

United Technologies Research Center large scale rotating rig and the first 
stage of the Space Shuttle main engine (SSME) high pressure fuel turbopump. 
Time-averaged blade midspan pressure and heat transfer profiles are calculated 
using the following different surface boundary conditions: adiabatic wall, 

prescribed wall temperature, and prescribed heat flux. Numerical solutions 
for the large scale rotating rig are compared with experimental data. Unsteady 
pressure envelopes are also presented for each geometry. In addition, instan- 
taneous contours are plotted for the SSME configuration. 
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1 . INTRODUCTION 


Current propulsion systems demand optimal performance and increased blade 
durability. To meet these demands, a clear understanding is needed of the 
unsteady aerodynamic and thermal environment associated with turbomachinery. 

In the past few years, turbine analysis capabilities have been enhanced through 
the development of computational fluid dynamics (CFD) codes that solve for the 
unsteady viscous flowfield in an axial turbine stage. The objective of this 
work is to extend and further evaluate the unsteady codes. Predictive capa- 
bility is assessed in terms of the accuracy of calculated aerodynamic and 
thermal blade loads and in terms of suitability for rocket propulsion 
applications. 

Specifically, R0T0R1, an unsteady rotor/stator interaction code, was 
modified to include heat transfer prediction capability. The original R0T0R1 
(an adiabatic wall version) and the modified R0T0R1 (a prescribed wall temper- 
ature or prescribed heat flux version) were applied to the midspan of two 
different turbine stages: the first stage of the United Technologies Research 

Center (UTRC) large scale rotating rig (LSRR) (for which tests were performed 
and reported in Refs. 1, 2, and 3) and the Space Shuttle main engine (SSME) 
high pressure fuel turbopump (HPFTP) turbine first stage. Results are pre- 
sented in the form of time-averaged blade pressure coefficients and heat 
transfer coefficients. All results obtained for the LSRR are compared with 
experimental data. In addition, time-varying flow features are also presented. 
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2. CONFIGURATIONS AND CONDITIONS 

The turbine stage flows analyzed in this study are described here. The 
first configuration considered is the single stage of the UTRC LSRR. This 
turbine stage has 22 stators and 28 rotors with an average blade axial chord 
of 15 cm (6 in.) and axial gap size of 15% average axial chord. Test con- 
ditions corresponding to the calculations performed in this study involved 
ambient air entering the turbine at 23 m/s (75 ft/s) and a mass flowrate of 18 
kg/s (40.3 lbm/s). The rotor speed was 410 rpm. The rig was run with and 
without a turbulence generating grid, and aerodynamic and heat transfer 
measurements were taken. Although the aerodynamics were not significantly 
affected by the presence of the grid, the heat transfer profiles measured with 
the grid in place were substantially different from those measured without the 
grid. Additional details are given in Refs. 1, 2, and 3. 

The second configuration analyzed is the first stage of the SSME HPFTP 

turbine. This turbine stage has 41 stators and 63 rotors with an average blade 

axial chord of 2.5 cm (1 in.) and an axial gap size of 0.8 cm (0.33 in.). The 

turbopump is driven by a gaseous hydrogen and steam mixture which, at full 

power level, enters the turbine at 226 m/s (741 ft/s), 70 kg/s (154 lbm/s), 

2 2 

1050 K (1900 R), and 37,500,000 N/m (5400 lbf/in ). The rotor speed is 
36,600 rpm. 
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3. METHOD OF SOLUTION 

This analysis was performed by applying an unsteady, viscous, single stage 
flow formulation to the midspan of the LSRR and HPFTP turbine first stages. 

The particular code used in this study is ROTOR1. 

3.1 CODE DESCRIPTION 

ROTOR 1 is the two-dimensional (2D) rotor/stator interaction code developed 
by M.M. Rai (Ref. 4). This code simulates the flow through an axial turbine 
stage by solving the 2D, unsteady, thin-layer Navier-Stokes equations. It 
features a factored, iterative, implicit algorithm, an Osher upwind differ- 
encing scheme, and a Baldwin-Lomax turbulence model. ROTOR1 employs multiple 
grids. Each blade is surrounded by a fine 0-grid which is overlaid onto a 
coarse H-grid. The H-grids are patched between blade rows, and the rotor 
H-grid slides past the stationary stator H-grid in time. ROTOR1 is a single- 
stage code and assumes that blade rows have an equal blade count. Reference 4 
provides a more detailed description of the R0T0R1 methodology. 

3.2 BOUNDARY CONDITIONS 

The inlet and exit boundary conditions used in this study are the same as 
those outlined in Ref. 4, so only a brief description will be given here. At 
the inlet, Riemann invariant R^ = u +[2/(Y-l)la is fixed and Riemann in- 
variant R = u -[2/(f-l)]a is extrapolated from the interior. The inlet 
flow angle is set to zero, and the entropy at the inlet is held constant at 
the freestream value. A constant static pressure is imposed at the exit. 
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Three different surface temperature conditions were explored: an adia- 

batic wall, a prescribed wall temperature, and a prescribed heat flux. The 
no-slip condition and a zero normal pressure gradient, along with one of the 
temperature conditions, comprise the surface boundary conditions. Reference 4 
describes the implementation of these boundary conditions assuming an adiabatic 

wall. If a nonadiabatic wall is assumed, the wall temperature, T , no longer 

w 

equals T^, the temperature at the first grid line off the wall. Implementa- 
tion of the surface boundary conditions assuming a prescribed wall temperature 
is described below. 

To impose no slip, the following equations are required: 

A (.p w u w ) = 0 
A (p w v w ) = 0 . 

To impose a zero pressure gradient, solve the energy equation 
K _ _e_ i ,(pu? 2 (pv? 2 , 

E - T-l + 2 1 p + p 1 


for P and take the normal derivative: 


_L_ _ ii _ 1 _9 f (pu) 2 (pv) 2 , 

Y-l 3n “ 3n 2 3n 1 p p J 


. 3p 

Setting 3 n = 0 and expressing the right-hand side in terms of discrete nodal 
points yields 


V V W,- V. 1 + W, - Vw> - 2 <V v w )( V p ,> • 

Applying the no-slip conditions and the time-change operator A gives 


2 <V V = ) 4p w + 4E W - 2 <V v w ,4 V “„ 4<p : , u a > + v „ 4 <p ,V- 4 V 0 
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To impose a constant blade temperature, use the perfect gas law and apply 



w 


T = 
w 


P T 
2 2 


Express T in terms of p , u , v , and E . 

2 2 2 2 2 


pT = [E-u(pu) 
2 2 c 2 2 2 2 

V 


v (P v ) + \ (u% v*)pj 

2 2 2 Z 2 2 2 


Apw - 


C T 
V W 


[( 


2 2 

U -4- V 
2 2 


)Ap - u A (p u ) - v A(p v ) + Ae ] = 
2 2 2 2 2 * 2 2 2 


The equations expressing no-slip, zero normal pressure gradient, and constant 
blade temperature can be written in matrix form. 


10 0 0 


A (p ) 

w 

0 10 0 


A (p u ) 
w w 

0 0 10 


A (p v ) 

w w 

a 0 0 1 

j 

A (E ) 
w 

— — | 

-B 6u 6v -6 

2 2 

. 

A (p ) 
2 

0 0 0 0 


A (p u ) 
2 2 

0 0 0 0 


A (p v ) 
2 2 

-a u V -1 

W W 


A (E ) 
2 


where 

a = (u 2 + v 2 )/2 

w w 

B = (u 2 + v 2 )/(2c T ) 

2 2 V W 

4 = l/(c v T w ) 
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The capability to calculate heat transfer was added to R0T0R1. The heat 
transfer was assumed to be due entirely to conduction, and therefore, Fourier's 
law of heat conduction was used, i.e., 



Two types of temperature distribution surrounding the blades were considered 

to approximate the normal derivative of T at the wall. First, the temperature 

at a blade surface point and at the two closest points along the normal line 

passing through the surface point was assumed to fit a quadratic polynomial of 
2 

the form a+bn+cn , where n is the normal distance from the surface. Next, 
the temperature distribution was assumed to be linear in n between the 
boundary point and the first point off the surface. Ideally, in both cases 
the off-wall point(s) through which the temperature polynomial is fit should 
lie in the viscous sublayer, i.e., y + should be no larger than 10. 

3.3 GRID 

For both the LSRR and the HPFTP, each 0-grid contains 151 circumferential 
grid lines and 51 radial lines. The stator H-grids have 66 grid lines in the 
l direction and 71 lines in the n direction, while the rotor H-grids con- 
tain 76 lines in the £ direction and 66 lines in the n direction. 

As mentioned previously, R0T0R1 assumes that blade rows have an equal 
blade count; in other words, there is a one to one ratio between stator and 
rotor blades. However, for the LSRR, the ratio between stator and rotor 
blades is three to four, and for the HPFTP, the ratio is two to three. In 
order to produce the appropriate mass flow rate in the calculation with just 
one stator and one rotor, the stator size must be reduced by 22/28 for the 
LSRR and 41/63 for the HPFTP, or the rotor size must be enlarged by 28/22 for 
the LSRR and 64/41 for the HPFTP. Two different grids were generated for the 
HPFTP, one with a reduced stator and one with an enlarged rotor. No dis- 
cernible differences in the solutions were found. The HPFTP grid using a 
reduced stator size is shown in Fig. 1. 
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4 . R0T0R1 ISSUES 


Several issues pertaining to ROTOR1 arose while performing the calcula- 
tions for this study. These issues are discussed in this section. 

4.1 REAL GAS VERSUS AIR PROPERTIES 


ROTOR1 assumes air to be the turbine's working fluid. This assumption is 
revealed in the use of gas properties in the code and the use of Sutherland's 
law to calculate p. The HPFTP turbine is driven by an l^/H^O mixture. The 
gas properties, such as specific heats, ratio of specific heats, and Prandtl 
number, can easily be changed from the air values to H^/H^O values. However, 
Sutherland’s law is not valid for the HPFTP gas mixture for any Sutherland's 
constant. To better simulate a non-air flow, a different way of determining 
the viscosity must be incorporated into R0T0R1. 

4.2 SURVIVING INITIAL TRANSIENTS 

When R0T0R1 begins its calculations using freestream conditions at all 
points in the computational domain as an initial approximation to the solution, 
strong transient pressure waves are generated. These initial waves can be 
strong enough to cause execution errors in the code. Several methods can be 
used to survive these transient waves. First, a small time step must be used 
at the beginning of the calculation. For the configurations in this study, 
R0T0R1 was run with 2000 iterations per blade pass for the first two blade 
passes (cycles). This number was decreased (thereby increasing the time step) 
to 200 iterations per blade pass by the seventh cycle. A way to lessen the 
strength of the initial transients is to bring the rotor up to speed slowly. 
R0T0R1 was run for 20 iterations with a rotor speed of zero. The rotor was 
then brought to full speed gradually over 500 iterations. If the transient 
waves still cause execution errors (as they did for the HPFTP configuration). 
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the eddy viscosity can be increased for several blade passes to damp the 
waves. The eddy viscosity was increased by a factor of 10 at the beginning of 
the calculation for the HPFTP . The viscosity was gradually decreased to its 
appropriate value over three blade passes. 

4 . 3 ITERATIVE PROCEDURE 

For some configurations and conditions, the solution to which R0T0R1 con- 
verges is inappropriate. Because a Riemann invariant, not actual values, is 
fixed at the inlet, the inlet values may change during the calculation. A solu- 
tion is then generated based on inlet conditions different from the desired 
conditions. In order to produce a solution with the correct inlet conditions, 
an iterative procedure is necessary. R0T0R1 is first run to convergence using 
the correct values of inlet Mach number and velocity. The converged value of 
inlet velocity is then compared with the correct value. If the velocities 
differ, the inlet Mach number is adjusted in the following manner: 

u u 

M . _ M + input - converged 

new old a . 

converged 

The new inlet Mach number, M , is then input into R0T0R1, becoming M , 

new old 

and the code is run until a periodic solution is obtained. Again, the con- 
verged inlet velocity is compared to the correct value, and the Mach number is 
adjusted if necessary. This procedure is repeated until the converged values 
of Mach number and velocity match the correct values. For the LSRR, the con- 
verged inlet Mach number and velocity matched the input values in the first 
run; no adjustments were needed. Calculations for the HPFTP were performed 
four times to obtain the correct inlet values. 
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5. RESULTS AND DISCUSSION 

Calculations were performed on the Cray X-MP 416 computer at Marshall 

Space Flight Center. A R0T0R1 calculation starting from freestream conditions 

6 

requires six CPU hours to reach periodicity and 2.1*10 words of core memory 
6 

plus 4.1*10 words of solid state storage device (SSD). 

Results generated using the modified R0T0R1 code are presented here. The 
results are divided into two classifications, aerodynamic and thermal. Aero- 
dynamic results are shown in the form of instantaneous contour plots, unsteady 
pressure envelopes, and blade pressure coefficients. The thermal results are 
presented in terms of heat transfer coefficients or Stanton numbers. 


5 . 1 AERODYNAMIC RESULTS 


Plots of calculated and experimentally measured pressure coefficients, 
defined as 

% ~ <P - P To> / ( 2 *0 " 2> 

are shown in Figs. 2 and 3 for the LSRR. These C^ profiles represent time- 
averaged values p, p ' , and p^. The agreement between the prediction and test 
data is excellent and consistent with results reported in Ref. 4. Pressure 
coefficients calculated for the SSME HPFTP are shown in Figs. 4 and 5. 

Figures 6 through 9 show the unsteady pressure envelopes for the LSRR and 
HPFTP. In both turbines, the effects of rotor/stator interaction, shown by 
the variations in pressure amplitudes, are experienced significantly by the 
rotor. Effects are also seen on the LSRR stator due to the narrow axial gap. 
The larger axial spacing in the HPFTP minimizes the effects of the blade 
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interaction communicated upstream in the turbine. Figures 10 through 12 show 
instantaneous contours for the HPFTP configuration. Figure 13 illustrates the 
flowfield unsteadiness which can be characterized by a ROTOR1 calculation. 

This plot shows the migration of the stator wake through the rotor passage. 

5.2 THERMAL RESULTS 


Stanton numbers for the LSRR were calculated using two different ex- 
pressions. One definition is consistent with the measured Stanton numbers 
reported in Refs. 2 and 3, i.e.. 


St = 


p U c (T - T_ ) 
e e p w To 


This definition poses a problem. As q" is not dependent upon T^ - T , but 

instead on T - T (Ref. 5), nondimensionalization by T - T yields a Stanton 
aw w To 

number dependent of T . The question of what wall temperature is the appro- 

w 

priate temperature to prescribe at the blade surface as a boundary condition 
becomes an important issue. Figure 14 shows calculated Stanton number for the 
LSRR rotor with three different wall temperatures. Three different Stanton 
number profiles result. 


The Stanton number was redefined to remove the wall temperature 
dependence : 


st = 


p U c 
e e p 


(T - T ) 
a w 


Figures 15 and 16 show blade-temperature independent Stanton numbers for the 
LSRR. Also shown are predictions resulting from specification of a constant 
heat flux, which are seen to be low in comparison to the prescribed wall temp- 
erature results. Experimental data are also shown. Comparison of the stator 
profiles shows the prediction to be essentially laminar on the pressure surface 
where the calculated profile shape resembles the grid-out measurements. Calcu- 
lations for the stator suction surface fail to correctly predict the transition 
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point or the profile shape. The calculated rotor pressure surface appears low 
and flat as compared with the experimental data. The shape of the suction 
surface profile, however, closely matches the shape of the grid in measure- 
ments, although the magnitude of the predicted profile is lower than that 
measured . 

Figures 17 and 18 show Stanton numbers for the HPFTP. Also shown on these 
plots are Stanton numbers calculated with a more conventional steady inviscid 
code, STAN5 (Ref. 6). A discrepancy between the R0T0R1 and STAN5 results is 
apparent. This could be due to the R0T0R1 assumption of air as the working 
fluid in its calculation of y. STAN5 incorporates real gas properties. 
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7 . CONCLUSIONS 


A numerical study of the unsteady aerodynamic and thermal environment 
associated with axial turbine stages was conducted, using the code R0T0R1. 
R0T0R1 was applied to two different turbine stages: the first stage of the 

UTRC LSRR and the SSME HPFTP turbine first stage. 

R0T0R1 is an effective tool with which to characterize the unsteady nature 
of flows through axial turbines. Instantaneous contours can be generated for 
any time within a blade pass, and unsteady flow features can be tracked through 
time. Time-averaged quantities can also be calculated with R0T0R1. Time- 
averaged pressure distributions calculated with R0T0R1 compare well with the 
experimental data. Qualitatively, the heat transfer profiles predicted for the 
LSRR exhibit reasonable agreement with the data, except on the stator suction 
surface where R0T0R1 fails to correctly predict the transition location and 
profile shape. The magnitude of the Stanton numbers is somewhat low, partic- 
ularly when a constant heat flux boundary condition is imposed in R0T0R1. 

The R0T0R1 -predicted heat transfer for the SSME HPFTP does not compare 
favorably with the conventional STAN5 prediction. The discrepancy is believed 
to be due to the ROTOR1 assumption of air as the working fluid for the vis- 
cosity calculation. This assumption yields erroneous heat transfer results 
for turbines that are not air driven. 
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Fig. 1 SSME HPFTP Turbine First Stage Grid 
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Fig. 2 Pressure Distribution on the LSRR First Stator 
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Fig. 15 Comparison of ROTOR1 Predicted and Experimental Heat 
Transfer for the LSRR First Stator 
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Fig. 16 Comparison of ROTOR1 Predicted and Experimental Heat 
Transfer for the LSRR Rotor 
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Fig. 17 Comparison of ROTOR1 and STAN5 Heat Transfer Predictions for 
SSME HPFTP First Stator (STAN5 Calculation by H. HcConnaughey 
of NASA-Harshall Space Flight Center) 
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Fig. 18 Comparison of R0T0R1 and STANS Heat Transfer Predictions for 
HPFTP First Rotor (STAN5 Calculation by H. McConnaughey of 
NASA-Marshall Space Flight Center) 
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